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ABSTRACT 


This paper evaluates five maneuvering strategies for an 
evader in a two dimensional continuous evasive game. An 
evader's movement, in two dimensions at constant speed, is 
Simulated by choosing courses that are independent and dis- 
tributed according to one of the five course change rules. 
The times between course changes are independent exponential 
random variables. An improvement over a previously estab- 
lished least upper bound on evader survivability is pre- 
sented. The bulk of this paper is the presentation of 
evader survivability as a result of the maneuvering strategy 


he employs. 
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I. INTRODUCTION 


A classic military problem with current application is: 
how to maneuver a mobile target in order to prevent success- 
ful prediction of its position. For example how should an 
aircraft carrier maneuver to maximize its survivability 
against an ICBM attack? 

Of primary importance to this problem is the time lag 
between the decision to fire at a target, made at time t, 
and the warhead detonation at time t+TL. This lag is called 
time-late and denoted as TL. The decision to initiate such 
an action is based in part upon knowing the target's posi- 
tion at time t while the effectiveness of that action is 


dependent upon the error in predicting that position at t+TL. 


A. ICBM VERSUS AIRCRAFT CARRIER 

Consider the ICMB versus aircraft carrier problem. As- 
sume that the carrier operates in mid-ocean, possesses great 
endurance, can make sharp turns and has as its only kinematic 
restriction a constant speed, v. The carrier is referred to 
as the evader, E. The attacker, P, continuously observes E 
from a nearby unarmed trawler. P's weapon is a land based 
ICBM which has perfect accuracy, produces a lethal area A 
(this lethal area can have any shape) but does not have a 


mid-course guidance nor homing capability. Although P knows 





E's position past and present, he also knows that his ICBM 
will have a one hour time-late.! 
P's problem is to predict E's position at time t+TL in 
order to initiate an attack at time t. E's problem, since 
he will not know of an attack until he observes the deto- 
nation, is to maneuver in such a way to confound P's predic- 


2 


tion. E's movement should be random because any nonrandom 


movement would be vulnerable to extrapolation. 


B. AN UPPER BOUND ON EVADER SURVIVABILITY 

A lower bound on P's success, Ры: and consequently an 
upper bound on E's survivability, l-ppr is known. P knows 
E's position at t+TL will be interior to or on a circle of 
radius v-TL centered on E's position at t (which is known). 
The enemy can then gaurantee a kill probability of at least 
A/1 (v: TL) © by choosing.a lethal area randomly within the 
circle of uncertainty (a wedge of random orientation will 
do) [Ref. 4]. Then the maximum survivability E could attain 
would be 1-A/m(v+TL) 2. 

This upper bound on E's survivability is a function only 
of P's strength, that is the magnitude of A (the value 
7 (v*TL)* is constant). Figure 1 is a graph of that survival 


function versus P's probability of success. 


trime-late is a summand of many factors in the command and 
control problem as well as the missle flight time. For the 
purpose used here it is sufficient to treat time-late in 
total and not consider its decomposition. 


“TF E knows the time at which an attack is initiated then 
the problem is trivial. 
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FIGURE 1 


Unfortunately the specific set of rules for E to maneuver 
mito achieve the survivability depicted in Fig. 1 is not 


known, 





II. BACKGROUND 


A. DISCRETE EVASIVE GAMES 

It was felt that a two dimensional continuous game was 
too difficult to solve. So the approach was to solve a 
similar but easier game. It was assumed that the ocean was 
a linear set of discrete points and that E's mobility con- 
sisted of being able to jump either right or left to an 
A jacent point [Refs, 2 and 3]. Time-late became an integer 
number of jumps that E could make between the attack de- 
cision and warhead detonation. As in the exact problem E 
had no knowledge of the decision to attack. 

Discrete problems of this type are called evasive games 
and are classified by the number of jumps constituting time- 
late. A one step discrete evasive game means E can move 
either left or right one jump prior to detonation. Game 
theory provides an immediate solution to that game; the 
value is p,=0.5. Also there is an optimal strategy for E 
in that game. At each jump E should go left with probability 
one-half or right with probability one-half. Employment of 
such a rule for each jump guarantees E a survival probability 
of at least the game value, regardless of P's firing rule. 

Having easily solved the one step game interest was 
focused on the two step game. In this game E was allowed 
two jumps during time-late. The solution has been obtained, 
but it was not as easy to achieve as that of the one step 


game [Ref.2]. 





The direction of the analysis was clear. Knowledge 
gained from solving the more simple games would be a step- 
ping stone to the solution of the more difficult games. 
Eventually the assumptions of linearity and discreteness 
could be relaxed. Progress, however, in solving the dis- 
crete games has been very slow. Researchers are presently 


embroiled in solving the three step game. 


в. A CONTINUOUS EVASIVE GAME 

Washburn in Ref. l presented a different approach to 
analyzing the ICBM versus aircraft carrier problem. He de- 
veloped the probability density function of a particle 
moving continuously in a two dimensional medium subject to 
a Specific set of maneuvering rules. That approach was 
distinct from previous work because it addressed the exact 
problem. The results presented in this paper are an ex- 
tension of that approach. 

Washburn proposed the following strategy for E and 
analyzed the evader's subsequent survivability, first as a 
function of the attacker's strength and secondly as a 
function of a specific strategy parameter. E was to select 
courses from a uniform distribution. Course changes were 
to occur as a Poisson process. This meant the time between 
course changes would be exponentially distributed. 

The exponential time between course changes has an 
intuitive appeal because of the memoryless property of that 
distribution. The memoryless property is: the probability 


that E will not change course in the future, given that he 
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has not changed course ЕЗ? some observation period, is in- 
dependent of the length of that period. Therefore, use of 
the exponential time distribution should serve to confound a 
prediction of future position by extrapolation, regardless 
of the course distribution used. For example, suppose P's 
attack decision rule required that he observe E maintain a 
constant course for at least five hours prior to initiating 
an attack. Such a procedure would not improve P's prob- 
ability of a kill because the only information of any benefit 
at the time an attack is initiated is E's position and last 
course. 

Figure 2 is a graph of the evader's maximum attainable 
Survivability when using the strategy of uniform courses 
and exponential times versus P's strength for the optimal 
exponential parameter, ЖС For comparison the upper bound 
on survivability is also presented. 

"If we assume the parameters A,v,TL are known to both 
sides, then the evader can select à to maximize the sur- 
vivability. The evader will clearly be in trouble if he 
makes Å too small, because the kill probability is at least 
exp (-I-TL).” On the other hand making course changes too 
frequently will lead to a density function that is highly 
peaked at the origin, which is equally undersirable..." 


[Ref. 1]. 


3The exponential parameter, A, is the inverse of the mean 
time between course changes. 


“Exp (-1+ TL) is the probability that E will make no course 
Ghange auring TL. 
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FIGURE 2 


"It is not known whether or not the optimal strategy for 
E is a Poisson strategy of the type just considered, or even 
whether the uniform distribution on angles is optimal within 


the class." [Ref. 2] 
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III. RESULTS OF SIMULATING THE CONTINUOUS EVASIVE GAME 


The goal of this research was to К ше if the uni- 
form course rule was optimal within the class of Poisson 
strategies or to find a better rule if it was not. For 
another rule to be better than the uniform one it would 
have to enable E to attain a higher survivability than that 
shown in Fig. 2. 

The methodology of solving discrete games had not ad- 
vanced sufficiently to achieve those goals. The methods 
used by Washburn were also not suitable because of the 
mathematical difficulty of the problem. 

The Poisson class of strategies is, however, uniquely 
Suited for analysis by computer simulation. This is due to 
the memoryless property of the exponential times between 
course changes. The results of such a simulation are pre- 
sented in this paper in the form of improved survivability 
for the evader. 

Five course change rules were evaluated in the simula- 
tion. The probability density function of each is pre- 
sented on subsequent pages along with a graph of the re- 
sulting attacker's py versus attacker strength, denoted as 
в-А/т(у-ть)2,7 These graphs are for representative values 
of the parameter a. Alpha, a, is the product of time-late 
and the exponential parameter i. 


Course changes were assumed to be independent of the 
underlying Poisson process. 
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The simulation Me a new lower bound on the p, 
the evader could yield to the attacker. Figure 3 is a graph 
illustrating the improvement of the new bound over that 
found by Washburn. This new bound is the least lower bound 
of all the p, versus $ curves for the four values of a and 
all rules simulated. Examples of such curves are shown in 
Figs. 5-9. The old bound in Fig. 3 is the lower bound of 
the py, curves simulated using the uniform rule exclusively. 

The minimal py, as a function of a, over the five rules, 
is shown in Fig. 4 for three specific attacker strength 
levels. The curves in Fig. 4 confirm the statement in Ref. 
2 that the evader should turn most often, that is a should 
be highest, against the weakest opponent. 

Comparisons of the Py Curves of thegdcifferente mies 
have shown that the uniform rule is not optimal for all 
attacker strengths and a values. For example, the reverse 
course rule (see Fig. 7) was shown to be better than the 
uniform rule for a=l. The difference between these two 
rules is illustrated in Fig. 10. Against a weak opponent 
(strength <.6) the evader would do better to use the re- 
verse course rule rather than the uniform rule for a=1. 
This improved survivability occurs because the reverse 
course rule created a "flatter" evader position density. 

None of the five course rules evaluated produced a con- 
Sistently smaller py for all values of attacker strengths 
and a. However, the left-right rule was dominated by the 


other four in all cases. Each of the five course rules was 
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simulated at six different levels of a by varying ET. The 
graphs of p, versus А/т(у»ті)2 for each simulation are pre- 
sented in the Computer Output section along with tabulated 
Pg S for specific attacker strength levels as a function 


of the course rule and the a used. 
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FIGURE 10 
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IV. SIMULATION MODEL 


Let E's position at time t be at X=0, Y=0 of a cartesian 
coordinate system. Then rotate the axis until the positive 
abscissa is aligned on and pointing in E's last direction 
of travel. Then regardless of the course rule used by E 
there is at least exp(-a) probability ENE he will be at 

=v-TL, Y=0 at time t+TL. This is the probability that E 
does not change course during TL. So the attacker's prob- 
ability of killing E by centering a small portion of his 
lethal area at (v:TL,0) is exp(-a). Then add to exp (-0) 
the integral of E's density function over that part of the 
uncertainty circle where the rest of the lethal area is 
targeted to determine the total p,. The total amount of 
lethal area allocated in this manner is the numerator, A, 
of the attacker's strength function S. 

A computer program to simulate the above computation 
was written in FORTRAN and run on the IBM 360. The program 
consisted of two phases, first the play of a strategy and 


secondly the scoring of that play. 


A. SIMULATION OF A STRATEGY 

To create a single play the evader's track was simulated 
from an initial position at time t to the resultant position 
at t+TL. The track was the result of a specific maneuvering 
strategy being simulated. A strategy was made up of two 
decision rules. The first rule determined the times be- 
tween course changes which were exponential random variables 


with mean ET=1/1. This rule was common throughout all of 


zu 





the simulations, although the parameter ET was a variable. 
The second rule determined the magnitude of subsequent 
course changes. The only other kinematic restriction was 
a constant speed, v, for the evader. 

To start a play datum was initialized by setting t=0, 
 X(0)=0, Y (0) -0 and an initial course, су, was selected from 
the uniform distribution. Then a sequence of exponential 


times [t;] were generated and a sequence of course 


6 
i=1,N (TL) 


changes, [er j L were generatéd using the course 


i=2,N(TL 


change rule. The evader's position at t=TL was 


N-1 

а) = D vet,» (Cos (c; )-Cos (cy) ) + v*TL*Cos (cy) IE 
N-1 

Y(TL) = ТЕ. vet, ° (Sin(c;)-Sin (cy) ) + veTL+Sin (cy) 1542 


В. SCORING THE PLAY 

Once E's position at t=TL was determined that observa- 
tion was scored. To score a play the coordinates of E's 
position, X(TL), Y(TL), had to be transformed. The purpose 
of the transformation was to make the observed position in- 
dependent of the particular initial cource су. The trans- 
formation was a rotation of the coordinate axis about the 
datum so the positive abscissa would be aligned in the 
direction of the initial course. The transformed position 


Was; 


othe number of elements of this sequence, N(TL), is a 
Poisson random variable. 
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XUL) 


N 


X (TL) Cos (c1) T Y (TL) Sin (c, I 2.1 


YTE) 


Y (TL)Cos (c,) = X(TL) Sin(c, ) 2„2 
An example of this transformation is shown in Fig. ll. 

If E had not made a course change during TL then his 
transformed position would have been 


X' (TL) 


l 
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A grid system of square cells’ was placed over the 


playing area of Fig. 11B and a determination was made as to 


А | 
I 
y X 
<A 
X 
RADIUS = v*TL 
ORIGIN: Datum ORIGIN: Datum 
Per), YTE) EU UNE). Gi) 
Ci: Initial Course A: (v*TL,O) 
SAMPLE PLAY PRIOR SAMPLE PLAY AFTER 
TO TRANSFORMATION THE TRANSFORMATION 
FIGURE 11A FIGURE LIB 


which cell X'(TL), Y'(TL) was in. Each cell of the grid 
had an associated value which would represent the number of 


times a play resulted in an observation in that cell. When 


The cumulative cell area was normalized by the factor 
T(veTL) 2, 
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the appropriate cell was determined, for the play being 
scored, that value was incremented by one. 

A simulation run was composed of 16810 plays and scoring 
iterations for a strategy utilizing a specific course change 
rule. The input variables for a run, besides the course 
rule were v, TL and ET. At the end of a simulation run the 
grid system was. a two dimensional histogram of E's position. 
The se in the histogram were then ordered, accumu- 
lated and normalized to achieve cumulative cell probabilities. 
The ordering corresponded to the conservative assumption that 
P could divide his lethal area and target only those cells 
with the higher probabilities. These probabilities were 
then plotted against the cumlative cell area they repre- 
sented. 8 This e was labeled p, versus attacker strength 
and was the primary output of the program. 

Thirty simulation runs were made to investigate five 
different course change rules and six different values of 
ET. For all runs the following constant values were main- 
tained, v-5 knots and TL-2 hours. The graphs from those 
runs are included in the Computer Output section. Each 
course change rule was simulated many different times, the 
only difference between runs being the pseudo random num- 
ber generator seeds. This was done to check for the vari- 


ЕЕ Бу of the p, graph for that rule. In all such runs 


The amount of area in a cell was denoted as cell size and 
egual to (2:TL/41)2. TL and v were held constant there- 
fore cell size was always approximately 1/4 nm“, 


24 


the resultant graphs were so similar that any difference 
was not distinguishable. For simulation run required two 


Minutes and forty-five seconds using 125 K on the IBM 360. 
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VV.  GONGEUS TON 


It was determined that the me Ec change rule 
is not optimal within the Poisson class of strategies. Com- 
parison of various rules, see Fig. 10, showed that under 
certain conditions a single rule such as the reverse course 
change rule is better. Also the new lower bound on p,, 
taken over all rules and values of a simulated, is an im- 
provement on the lower bound achieved from only the uniform 
mule. 

Certainly not all the possible course change rules were 
Simulated. The rules evaluated, however, were representa- 
tive of the broad class of possible rules. When compared 
to the uniform rule, all but one of the other rules showed 
that E could improve his situation if he knew P's strength 
by selecting the better rule for that encounter. The one 
rule that was consistently dominated was the left-right 
rule. 

It remains unknown whether or not the optimal strategy 


1s a Poisson strategy. 
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APPENDIX A 


FLOW DIAGRAM 


START 


IX -- 55555 --- ipseudo-random number seed, an input variable | 
= Evade: speci DLnots | 
IL == 2.0 === time late Mour S | 


ШІ -- 1.0 --- mean time between course changes, an input 
"variable 

BE(41,41) --- playing area and histogram; each cell is 

PP(i,j) -- 0 la square with dimension (2*V*TL)/41 


mor all i,j 
---ı INDEX is a book-keeping record of PP cell 
° ‘feasibility 






INDEX (41,41) 
INDEX(i,j) -- 0 
mor all i,j 


INDEX(i,j) -- O see PE(i j) an exterior сеш 
| No 
INDEX(i,j) -- 1122415 PP(i,j) an interior cell 


No 
INDEX 2 


© 





х -- 0.0 
y =- =0n0 
Clock -- 0.0 
c -- 0.0 


TR -- TL-Clock (3) 
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E Generate c=U(0,2 ) 
Theta -- c 


Generate twExp(1/ET) 
T -- min(TR,t) 


Clock -- Clock +t 


- 








Generate a course 
change 












x -- x + (V4TL)*Cos(c) 
y -- y + (V*TL)*Sin(c) 


| TR == Th=Clock No Is Clock = TL 


Yes 






х -- x*Cos(Theta)t+y*Sin (Theta) 
y -- y*Cos(Theta)-x*Sin(Theta) 











Determine which column of PP the x co- 
ordinate is in, e.g. column 1 








Determine which row of PP the y coordi- 
nate is in, e.g. row m 


PP(m,1) -- PP(m,1) + 1 






Sort and order PP 


Sum PP 


Graph: p, versus attacker strength, 
STOP М k À 
_—— i.e. py versus normalized lethal area 





Subroutine Trktm(ET,r,t) 
t -- -ET*ln(r) 







Subroutine Corse(r,c) 


igen io. Lal) 


DP 


delta c -- r*(b-a) + c 


C -— c t delta c 


en 
| 
! 
Fh 


с --- с,Моа(Ғ) 


Return 


Main 


Subroutine Tcorse(r,c) 


s -- 0.2 
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No 


Pi -- 3.16 


т=— жир 


Is r between О апа 0.25 inclusive 


Yes 


a -- s/2 


b -- 1/(2*Pi) - s*Pi/4 


с. 
| 
| 
| 
tok 
Q 
> 


e -- (b**2)-4*akd 


delta CZ NNN 


Is r greater than 0.25 and less than or 
equal to 0.5 


Yes 


a -- -s/2 





Рр Same Pi) Зехра 


а -- (-s*(Pi**2)/4) - r 


eu GEN E 


delta c -- (-b + Sqrt(e))/(2*a) 





No 


Is r greater than 0.5 and less than or 
equal to 0.75 


Yes 


rprime -- r - 0.5 
a -- s/2 
b -- 1/(2*Pi) - s*Pi/4 


d -- -rprime 


е -- (b**2) - 4%a%d 
delta c -- (-b = Sart(e))/(2*a) 


delta c -- delta c + Pi 


rprime == Lan 


b -- 1/(2*Pi) + 3*s*Pi/4 


d -- (-s*(Pi**2)/4) - rprime 


ш 

| 

! 

N 
Qo 
~~ 
Мә 


е —— (Ь**2) - 4*akd 
delta c -- (-b * Sqrt(e))/ (2*a) 


delta c -- delta c + Pi 


(А) 
H 








c -- € + delta c 


c -- c,Mod(f) 






Subroutine Pcorse(r,c) 


РІ -- ЗГЕ» 


f -- 2*Pi 


s -- 1/(Pi**2) 


No 
Is r greater than 0-53 
a -- -s/2 delta c -- Sqrt(a) 


b -- 1/(2*Pi) + 3*Pi*s/2 


| 


| | 
(v 
o 


d -- -s*(Pi**2) - r 


e -- (b**2) - 4*akd 


delta cy -- (-btSqrt(e))/ (2*a) 





delta е7 а (-b-Sart(e))/(2*a) 


D 





Is Pi less than delta cj and is 
delta cı less than or equal to f | Yes 


delta c -- delta с» delta c -- delta сі 


| 
О 






c-- c+ delta c 








c -- c,Mod(f) 


Subroutine Treners(r,c) 


Pi ONS 


f -- 2*Pi 
а -- Р1/6 
b -- 11*Pi/6 


delta c -- r*(b-a) +a 


c -- c+ deltac 


c -- c,Mod(f) 


с. 
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C -- c + Pi/2 


No 









Pi -- 3.1415 








Yes 


Subroutine Lrcors(r,c) 


Is r greater than 0.5 





Rule 


Uniform 
Modified 
Left-Right 


Reverse 
Course 


Left-Right 
Truncated 
Uniform 


column 
minimum 


Rule 


Uniform 
Modified 
Reverse 
Course 


Left-Right 


Truncated 
Uniform 


COMPUTER OUTPUT 


TIET 
10 2-0 3.0 4.0 
129 . 280 ‚501 ‚486 
797 «580 2501 ‚486 
‚690 .572 ‚548 ‚564 
917 ‚768 ‚650 ‚588 
‚ 697 2250 .309 25 1 
. 690 ‚556 ӘЛЕ .486 
Py for attacker strength equal 0. 
Table I 
TL/ET 
ДЕЛІП 20 3.0 4.0 
‚854 ‚784 ЛЭХ ‚760 
‚878 . 784 2752 . 768 
Left-Risht 
.838 . 791 „799 ТОВ 
.964 .909 „870 „854 
‚846 . 768 . 768 815 
.838 . 768 122 ‚760 


column 
minimum 


рр for attacker strength equal 0.4 


Table II 
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TL ET 





Rule 1.0 20 350 4.0 
Uniform ‚933 . 909 . 901 „917 
Modified ‚948 . 909 . 901 .917 
Left-Right 

Reverse . 933 .917 .925 . 948 
Course 

Left-Right «987 «987 ‚980 „987 
Truncated ‚933 . 909 . 909 . 933 
Uniform 

column 933 .909 . 901 . 917 
minimum 


Pk for attacker strength equal 0.6 


Table III 
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